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Abstract 

Making use of the Effective Field Theory(EFT) expansion recently developed by the authors, we compute 
the charge form factor of triton up to next-to-next-to-leading order (N^LO). The three-nucleon forces(3NF) is 
. required for renormalization of the three-nucleon system and it effects are predicted for process and is quali- 
^ tatively supported by available experimental data. We also show that, by including higher order corrections, 
the calculated charge form factor and charge radius of are in satisfactory agreement with the experimental 
data and the realistic Argonne vis two-nucleon and Urbana IX potential models calculations. This method 
f^piakes possible a high precision few-body calculations in nuclear physics. Our result converges order by order 
^ in low energy expansion and also cut-off independent. 
■ PACS numbers: 21.45. -Fv, 25.10.-Fs, 25.20.-x, 27.10.-hh 
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I. INTRODUCTION 



The few-body problem has now reached the stage when the main interest is centered not on theoret- 
ical techniques but on the application of these techniques to problems in nuclear, particle and atomic 
physics. Among these theories, EFT is a powerful tool to calculate low-energy observables in a system- 
atic way. Pionless EFT has been applied to two- and three-nucleon systems H, 0, H, 0, 0, Isl, 0] , 
during the last few years. It is ideally suited to exploit a separation between scales and can use with 
contact interactions at low-energies. Some of these low-energy observables such as neuron-deuteron 
radiative capture and the two-body photodisintegration of the triton have been calculated using pio- 
nless EFT and insertion of the three-body force, up to N^LO [l^, 11, 12, 1^. In these calculation we 
find a simple way for inclusion of external currents. The evaluated cross sections have been compared 
with experimental data and other modern realistic two- and three-nucleon forces AV18/UrbIX poten- 
tial models calculations. We are also going to study the structural functions of three-body systems 
in future calculations. Among them, the study of the three-body nuclear physics involving the charge 
form factors and charge radius of is calculated by using pionless EFT. It have been investigated in 
theoretical and experimental works over the past decays. 

The effect of meson exchange currents on the charge and magnetic form factors of was investigated 
by Maize et al. [13] and ref there in. Inclusion of meson exchange currents considerably improves the 
impulse approximation fits to the experimental data. An other calculation reported in include single-A 



isobar admixtures in the three-nucleon wave function [15 

Recently, a very good agreement reached with the calculation of the charge form factor of '^H 16|, |l7 . 
They calculated the electromagnetic form factors of the three-nucleons with wave functions obtained 
with the Argonne vl8 two-nucleon and Urbana IX three-nucleon interactions. The two-body currents 
required by current conservation with the vl8 interaction as well as those associated with NA transition 
currents and the currents of A resonance components in the wave functions, are considered for these 
calculations. In this way, explicit three-nucleon current operators associated with the two-pion exchange 
three-nucleon interaction arising from irreducible S-wave pion-nucleon scattering is constructed and it 
is shown to have very little effect on the calculated magnetic form factors. 

More recent calculation of the triton charge form factor at leading order in (i/\a\ and extract the 
charge radius has been performed by Platter et al. 1^ . They studied the correlation between the triton 



binding energy and charge radius and relate it to the three-body parameter L3 that parameterizes the 
Phillips and Tjon lines in effective theory. 

In experimental works, the charge form factor was measured in the range 0.29 < < 1 fm~^ 
by Beck et al. jioj. After few years, the charge and magnetic form factors were also measured in the 



range 0.0477 < < 2.96 fm ^ 20|. The experimental result of triton charge radius is quoted as rc= 
1.63 ± 0.03 fm [20]; as = 1.76 ± 0.04 fm Q; as = 1.81 ± 0.05 fm [23(; and as = 1.755 ± 0.086 
fm Q. 

In this paper, we study the charge form factor and chargr radius of the triton to next-to-next-to- 
leading order (N^LO). Our results show that, by including higher order corrections, the calculated 
charge form factor of is in satisfactory agreement with the modern nucleon-nucleon potential AVI 8 
together with the three- nucleon force UrbanalX. This result has several implications. The first and 
most important one is the remarkable success of EFT improving for simple inclusion of external currents. 
The second is the indication that the model makes possible a high precision three-body calculations in 
nuclear physics. It now appears to explain the structural functions of three-body systems. The third 
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is remarkable quality of the EFT low-energiy calculations. The result converges order by order in low 
energy expansion, model independent and also cut-off independent. It should be considered that it is 
a quantitative question based on current choices and gauge invariance of nuclear force models to reveal 
signatures by switching on and off 3N forces. 

This paper is divided into five main sections. Section II contains a brief description of the calculation 
of the the bound three-body systems. Faddeev integral equation, Lagrangian and kernel for the reaction 
involved in our consideration and three-body forces are given in this section. In Section III, we calculate 
the charge form factor of the triton by using EFT. Calculation of the charge form factors, the relevant 
diagrams, and parameters are discussed in this section. In section IV, comparisons of our results with 
the corresponding experimental and theoretical data are given. Finally section V contains a concluding 
discussion. 




FIG. 1: Faddeev equation for the three-body amplitude including the three-body force. Line shows nucleon 
and thick solid line is propagator of the two intermediate auxiliary fields dg and dt- The three-body interactions 
with strength Ti, denoted by the dot. 



II. BRIEF REVIEW OF THE BOUND THREE-BODY SYSTEM 

The derivation of the integral equation describing neutron-deuteron scattering has been widely 
discussed jHI, [9|. We present here only the result, including the new term generated by the two- derivative 
three-body force. The three- nucleon lagrangian is given by [6|: 

C = (ido + Si)n + df {-ido -Sl + ^s)df + (-i9o - S + A,) dl (1) 
+tt (i9o + g + 2! + t _ (rff (iV^P^iV) + H.c.) - g, (rff (iV^PW) + H.c.) 
-ujs {t\T^N)df + H.c.) - ujt {t\a'N)4 + H.c.) + . . . , 

where is the nucleon iso-doublet and the auxiliary fields t, d^ and d\ carry the quantum numbers 
of the '^H spin and isospin doublet, ^5*0 di-nucleon and the deuteron, respectively. The projectors are 
P* = -^T20'2(J^ and = -^(72^2^^, where A = 1,2,3 and i = 1,2,3 are iso-triplet and vector indices 

and (o"') are isospin (spin) Pauli matrices. 

For calculation in this channel, two amplitudes get mixed(see fig. 1): tg describes the dt + N — > dg+N 
process, and U describes the dt + N ^ dt + N process: 
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with 



tsip,k) =l[3}C{p,k) + 2n{E,A)] + ^ jdqq' [V^) [K{va) + '2^'H{E , K)]Uq) 



+Vt{q) [3/C(p,g) + 27^(E,A)]t,(g)] 

A 

U{p, k) = i [/C(p, k) + 2niE, A)] + lJ dq q' [D,{q) [/C(p, q) + 2n{E, A)] t,(g) 



+V,{q) [3/C(p, q) + 2n{E, A)] , (2) 



-Dsiq) = ^ + (3) 



4^ 



-7.+\/i(92-fc2)+7? 1 I -7.+;/i{g2-fc2)+7? 
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where Vs^tiq) = T^s,t{E — 9) are the propagators of deuteron. 

For the spin-triplet S-wave channel, one replaces the two boson binding momentum 7 and effective 
range p by the deuteron binding momentum 7^ = 45.7025 MeV and effective range pt = 1.764 fm. 
Because there is no real bound state in the spin singlet channel of the two-nucleon system, its free 
parameters are better determined by the scattering length = l/7s = —23.714 fm and the effective 
range Vg = 2.73 fm at zero momentum. The neutron-deuteron J = 1/2 phase shifts __6 is determined by 
the on-shell amplitude tt{k, k), multiplied with the wave function renormalization 



Tik) = Ztt{k,k). (4) 

In that case, a unique solution exists in the ^S'i/2-channel for each A and vanishing three-body force, 
but no unique limit as A ^ 00. As long distance, phenomena must however be insensitive to details 
of the short-distance physics (and in particular of the regulator chosen), Bedaque et al. [g], 0] showed 
that the system must be stabilized by introducing of a three-body force in a form: 

00 

m A) = A ^ ij,„(A) (H£^2i)^ (5) 

n=0 

which absorbs all dependence on the cutoff as A — > cxd [6.]. At leading order, 7 <^ A, the essential 
observation are of order 0{Q/A) and are independent of energy (or momentum) and hence can be 
made to vanish by only Hq and not any of the higher derivative three-body forces. At NLO, there is 
only a perturbative change from the LO asymptotic of order (!?(1/A^)(see |^). 

For N^LO order, a new three-body force seems to be required. The terms will be proportional to 
ME arising from expanding the kernel in powers of q/A and V ME/ A and a three-body force term 
which contains a dependence on the external momenta k"^ and ME can absorb them. It means we 
need H2 for calculations up to N^LO. It has been shown that the three-body system at N^LO can be 
renormalized without the need for an energy dependent three-body force at this order |8|. 
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In the presence of an electromagnetic interactions, after performing minimal substitution, requires 
that V operators be replaced by covariant derivatives D = V — ieA, and time-derivatives do be replaced 
by Dq = do + ieAo in lagrangian eq. (1), where A is the vector potential. 

This simple perturbative expansion of the scattering amplitude is reproduced order by order in the 
pionless expansion. Relativistic corrections are encountered at N^LO and effective range expansion 
provides a complete description of scattering in the low-energy region]^. 




FIG. 2: Diagrams that contribute to the triton charge form factor. Wavy line shows photon. Remaining 
notation as in fig. 1. 



III. THE CHARGE FORM FACTOR OF TRITON 

The electric charge form factor of the triton is measured over a wide range of momentum transfers. 
The modern nucleon-nucleon potential AV18 together with the three-nucleon force UrbanalX potential 
models and effective theory reproduce the data very well in the kinematic regions where they are 
applicable. (For a recent review, see [13, fis'].). In present work, we calculate the charge form factor in 
the range < < 0.3 fm~^ in corresponding to recent effective theory calculation 



A convenient expression to calculate the charge form factors of a J=l/2 nucleus, such as the triton, 
is obtained by the matrix element of the electric current, is given by: 

{^+\Je\^+) = eFciq') , (6) 

where q is the momentum transfer and is the normalized ground state wave function. The dimen- 
sionless charge form factor defined in eq. ([6]) are normalized such that Fc{0) = 1. The slope of the 
charge form factor at low-momentum transfer defines the triton charge radius: 

W) = i-q'(^')/6 + .... (7) 

where rc = (r^)^/^ is the charge radius. 

In the pionless EFT, the charge form factor Fc{q'^) has an expansion in powers of Q, Fc{q^) = 
F^\q'^) + F(^\q^) + FQ\q^) + where the superscript denotes the order of the contribution. 

The diagrams that contribute to the deuteron charge form factor in pionless EFT are shown in 
fig. 2. In addition to diagrams where the photon couples to the nucleon, there are also couplings to the 
dibarion field obtained by gauging the Lagrange density in eq. (1). 

The LO order calculation is corresponding to the leading term in a effective range expansion and 
is in corresponding with the charge form factor computed in ref. ll]. At higher orders there are 



contributions from higher dimension operators involving more derivatives on the nucleon field, such as 
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the nucleoli charge radius operator, and also from higher dimension couplings of the dibaryon field, 
that is the dibaryon charge radius. At NLO, we have contributions coming only from insertions of the 
kinetic energy in the deuteron propagator. The NLO corrections to neutron-deuteron scattering are 
simply the range corrections in first order perturbation theory. 

At N^LO order there will be corrections to neutron-deuteron scattering which are simply two times 
insertions of the kinetic energy in the deuteron propagator as well as correction from nucleon structure 
comes from the nucleon radii [2I, [^. For higher orders in the EFT expansion, order Q^, there are 
contributions from several types of graphs, including relativistic corrections and contribution from, a 
four-nucleon-one-photon operator, the dibarion charge radius operator. 

The last diagram also show the photon couples to the three-nucleon force. At LO calculation, there 
is no interaction of Hq with a photon, because Hq has no derivatives, so it is not affected by the minimal 
substitution P ^ P — eA. In other hand, contribution from the photon couples to the three-nucleon 
field (H2) is calculated at N^LO. 

Low-energy observables must be independent of an arbitrary regulator A up to the order of the 
expansion. One can therefore estimate sensitivity of the results to cut-off A, and hence provide a 
reasonable error-analysis, by employing a momentum cut-off in the solution of the Faddeev equation 
and varying it between the breakdown- scale A^ and 00. The cut-off varied between A = 200 MeV and 
A = 550 MeV. 



IV. RESULTS AND DISCUSSION 

A numerical solution of the charge form factor and charge radius of triton have been done. The 



numerical values and parameters used for calculation is described in ref. We solved integral 

equation by insertion of Q in Faddeev equation where is shown in figs. 1 and 2. The integral equation 

is solved by the first two terms of the effective range expansion, properly iterating and folding to 

electromagnetic interaction where are shown in fig. 2 order by order and integrated on the involving 

momentum. The three-body parameter Hq is fixed from the ^Si scattering length 03 = (0.65±0.04) fm 

2 

at LO. H2 is also required at N^LO and it is determined by the triton binding energy = 8.48 MeV. 

The calculated charge form factor is compared with the experimental data. The results have 
been shown in fig. 3. The rhomboids and squares are experimental data from Refs. [2l| and [l9|, respec- 
tively. While the general trend of the low-momentum transfer form factor data is reproduced by our 
calculation, we obtain a somewhat close slope to the experiments up to N^LO. The remarkable success 
of the present calculation by using EFT should be stressed. The number of three-body calculations 
with external currents is however, extremely limited. It suggests, in particular, that the present model 
provide simple inclusion of external currents. Yet, the excellent agreement between the calculated and 
measured results of charge form factor suggests that these corrections may be negligible. 

Finally, values for the charge radii of are listed in Table I. The results in comparison of the 
modern nucleon- nucleon potential AVI 8 together with the three- nucleon force UrbanalX, are found to 
be in good agreement with experimental data. Fitting a polynomial in to our result for the charge 
form factor, eq.(7), we can extract the triton charge radius. This procedure leads to a triton charge 
radius =[1.959(LO)-0.113(NLO)-0.071(N2LO)]fm =[1.775 ± 0.02] fm. 

The cut-off variation has been also done. These variations were nearly independent of variation of 
momentum and decreased steadily when the order of calculation is increased up to N^LO. 
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FIG. 3: The triton charge form factor. Dot, dash and sohd curve are LO, NLO, and N^LO EFT results. The 
rhomboids and squares are experimental data from Refs. |20i] and [2l|, respectively. 



TABLE I: Comparison between different theoretical and experimental results for triton charge radius in fm. 



Theory 


Triton charge radius 


% Error 




Effective theory [18] 


2.1 


20 


Arj 


ronne vl8+Urbana 


IX [17] 


1.722 


2 


Arj 


ronne vl8+Urbana 


IX [16] 


1.725 


2 




EFT(LO) 




1.959 


12 




EFT(NLO) 




1.846 


5 




EFT(N2L0) 




1.775 


1.1 




Exp. [24] 




1.755±0.086 





In fig. 4, we show the correlation between the triton charge radius and binding energy of triton. The 
solid line denotes leading order effective theory result 18|. The two boxes show modern calculations 



(including meson exchange currents) based on the AV18 potential with and without the Urbana IX 

The circles indicate Faddeev calculations using different potentials from Refs. 

The triangle with error-bar 



three-body force [17 

The correlation curve of Ref. 
indicates the experimental values j23 ]. 



27j is given by the dotted line. 



Short-dash, dot-dash and long-dash curves are LO, NLO, and 



N LO our EFT results, respectively. It should be considered that the inclusion of the range correction 
moves the universal line closer to the experimental values in corresponding to Phillips line. 
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FIG. 4: The correlation between the triton charge radius and binding energy. The solid line denotes leading 
order effective theory result 18(]. The circles indicate Faddeev calculations using different potentials from Refs. 



251. 126 1 w hile the two boxes show modern calculations using AV18 with and without the Urbana IX three-body 



force 17 1 . The correlation curve of Ref. [27] is given by the dotted line. The triangle with error-bar indicates 
the experimental values 2J]. Short-dash, dot-dash and long-dash curves are LO, NLO, and N^LO our EFT 



results, respectively. For better comparison, axes scales are chosen like as 18l |. 



V. SUMMERY AND CONCLUSION 

The goal of this work was to obtain the charge form factor and charge radius of triton, as a three-body 
bound state, that are obtained by using EFT. We applied the formalism developed to the computation of 
processes involving external currents such as, triton charge form factor and charge radii. We considered 
the diagrams needed to be computed and the charge form factor numerically obtained by perturbatively 
solving of Faddeev equation, up to N^LO. At very low energies, the interactions between nucleons can 
be only described by point-like interactions. The present study of the triton charge form factor show 
that it is possible to obtain remarkably good predictions for electrogenic structure functions of the 
few-body bound states by using EFT at very low-energies. 

The charge form factor of triton at low energies was calculated by using pionless EFT and simple 
inclusion of external currents. The charge radii is also calculated up to N^LO. The triton binding 
energy and nd scattering length in the triton channel have been used to fix them. Hence the triton 
charge radii is in total determined as =[1.959(LO)-0.113(NLO)-0.071(N2LO)]fm =[1.775 ± 0.02] 
fm. It converges order by order in low energy expansion and is also cut-off independent at this order. 
The results also move the universal line closer to the experimental values when order of calculation 
increases. This is in corresponding to the Phillips line. 
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